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Preliminary limits on the D° mixing parameter y = Ar/2r are obtained using about 57.8 fb 1 
of data collected by BABAR in 2000 and 2001: y = (1.4±1.0(stat.)lo.?(syst.))%. y is extracted, 
provided that CP is conserved, by measuring separately the D° lifetime for the Cabibbo- 
suppressed decay modes K~ K + , tt~-k + and the Cabibbo-favored mode K~tv + . Backgrounds 
are suppressed by D*-tag and particle identification requirements. 



1 Introduction 

Mixing in the charm sector is a suitable place where effects due to physics beyond the standard 
model (SM) can show up. The D° — D° mixing parameters are defined as x = AM /T and 
y = AT/2T with AM = M 1 - M 2 , AF = T 1 - T 2 and T = Ti + T 2 /2, and where M 1)2 an d 
T\ : 2 are the masses and the decay widths for the physical eigenstates |-Di,2) = p\D°) ± q\D°). 
In the SM these parameters are predicted to be very small and, even including all the proper 
corrections, one gets \x\, \y\ < 1CP 3 . However non-SM processes can significantly enhance 
x whereas final state interactions and SU(3)-breaking can enhance y to a level accessible by 
the current experimental sensitivity (10 -2 ) or in the near future at asymmetric B-factories 
(few x 10~ 3 ). 

One possible search for mixing effects can be performed by measuring the lifetime difference 
between D° decaying, through a Cabibbo-suppressed diagram, to a CP-even eigenstate (such 
as K~ K + and 7r~7r + ) and decaying to the Cabibbo-favoured K^tt + final staten. Assuming the 
latter to be an equal mixture of CP-even and CP-odd eigenstates the rate asymmetry for neutral 
D decays into CP + and CP~ eigenstates can be expressed in terms of a lifetime ratio: 



fjCP+) - f(CP-) _ t(K-tt+) 
f(CP+) +f(CP-) ~ r(h-h+) 



1 , h~h + = K~K + ,-k-^ 



where r = 1/T and T is the effective decay rate obtained by fitting to a pure exponential the 
time-dependent rate for each final state. The parameter y CP can be expressed!, by a good ap- 
proximation, in terms of the CP-violating weak phase (ft and \q/p\- Ucp ~ycos(f)—l/2\q/p\xsm(j). 
Thus, in the CP-conserving limit, y = y CP . 

Because of the similar topology of the above final states, many systematic uncertainties in 
the D° lifetimes cancel in the ratio so that the lifetime ratio provides a particularly sensitive 
measurement of y. 



2 Data sample and the BaBar detector 

The preliminary measurement of y presented here is based on a data sample of 57.8 fb _1 collected 
with the BaBar detector, at the PEP-II asymmetric e + e~ collider, during the years 2000 and 
2001. This sample includes data taken both on and off the T(4S) resonance since the analysis 
selects events from the cc continuum. This analysis uses also Geant4 simulated data, either 
generic qq samples of about 30 fb -1 and signal samples of various luminosities. 

Detailed description of PEP-II storage ring and the the BaBar magnetic spectrometer can 
be found elsewhere 13. The BaBar subdetectors employed in this analysis are the silicon vertex 
tracker (SVT), the drift chamber (DCH) and the ring- imaging Cherenkov detector (DIRC) that 
provide vertexing, tracking and K/tt separation. 



3 D° candidates selection and proper time reconstruction method 

D° candidates were selected by looking for pairs of charged tracks (h + h~ ) with combined invari- 
ant mass close to the D° mass. Each of them was required to satisfy track quality criteria for a 
good reconstruction. The x 2 probability of the common vertex to which the h~ were fitted 
was required to be better than 1%. Particle identification criteria were applied to both the tracks 
to distinguish among the three decay modes and to suppress background. The kaon selection 
was rather tight and characterized by a pion contamination of less than 3% for momenta less 
than 3GeV/c. The pion selection included a muon veto. Additional rejection of combinatorial 
background due to low momenta pions was obtained by cutting on the angle of a pion track in 
the D° center-of-mass with respect to the D° flight direction. 

The D° candidates were required to be produced by D* + — > D°irf decay. This allowed 
the introduction of the difference between the D* and the D° reconstructed mass, 5m = 
m(h + h~ir s ) — m(h + h~), which is crucial for background rejection. To increase acceptance, 
ir s candidates were required not to contain DCH hits but to contain at least 6 SVT hits. 

The selection of D* candidates coming from cc continuum, namely the rejection of those 
produced in B meson decays, was obtained by asking for a D° momentum in the e + e~ center- 
of-mass frame greater than 2.5GeV/c. This constrains the D* to be within the interaction 
region (beam spot) thus allowing a refitting technique with beam spot constraint which highly 
improves 5m resolution. The D* decay point was located by pointing the D° momentum vector 
back to the beam spot and the slow pion, which suffers from multiple scattering, was then 
refitted to this point and the D* itself obtained by a vertex fit involving the refitted ir s . The x 2 
probability of this fit was required to be better than 1%. The <5m-cut was performed by opening 
a ±2 -^3MeV/c 2 window around the 5m distribution (Figure []]) peak, depending on the quality 
of the 7r s track. Figure |l| shows the final invariant mass distributions; there were found about 
158, 000 D° -» Ktt, 16, 500 D° -> KK, and 8, 350 D° -> vrvr candidates with signal purities of 
about 99.5%, 97.1% and 92.4% respectively, within the D° signal region (indicated by dashed 
lines in Figure |l|) . 

The D° proper time is derived, in the plane transverse to the beams direction, from the flight 
length defined as the distance between the D° and D* vertices projected on the D° transverse 
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Figure 1: On the left the reconstructed D invariant mass distributions for the three decay modes after the 
application of the selection criteria. Satellite peaks are due to the reflections of D° — » K~n + decays deriving 
from incorrect mass assignment. On the right the corresponding 8m distributions concern, apart from the <5m-cut 
itself, the selected D° candidates falling within the mass window indicated by dashed lines into left plots. 



momentum. The choice to measure the flight length in two dimensions derives from the particular 
beam spot transverse flatness (o~ y ~ 5fim, a x ~ 120/im, o~ z ~ 9000/um) in connection with the 
use of the beam spot constraint. 



4 D lifetime extraction 

An unbinned maximum likelihood fit was used to extract the lifetime for the three D° samples. 
The likelihood function is divided into two different decay time distributions, one for the signal 
and one for the background. The former is a pure exponential whereas the latter is the sum 
of an exponential for the flying (charmed) background and a 5 function with a free offset for 
the not-flying (light quark composed) background. Each distribution is smeared by convolution 
with a proper time resolution function. 

The resolution model, successfully tested for simulated data samples, is a sum of three 
gaussians in which the first two have a width proportional to the event-by-event proper time 
error derived from the vertex error matrix and the third, of fixed width, is intended to describe 
tail contributions. A fourth wide gaussian is added in the smearing of the background function 
in order to account for long tails. 

To combine the signal and background likelihood functions, the D° reconstructed mass 
distribution was separately fit to derive the probability for each D° selected candidate to belong 
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Figure 2: The white histogram represents the result of the unbinned maximum lilelihood fit described in the text; 
the gray one is the portion assigned to the background by the fit. 

to the signal. The sidebands were included in the fit region to better constrain the background. 

The results of the lifetime fits are shown in Fugure ^ superimposed on the proper time 
distributions for the three decay modes. 

5 Lifetime ratio systematics 

D° lifetime measurement is a high precision measurement since the statistical uncertainty for 
the Cabibbo-favoured sample is about 1.3 fs (3 parts per thousand). At this level of precision 
wide studies of systematic sources (tracking, vertexing, alignment) are needed. 

At this stage we do not present an absolute D lifetime value because the overall systematic 
uncertainty (about 3-=-4fs) can be reduced with further studies. We checked the full consistency 
of the fit result with both our previous check § on a sample of 12.8 fb -1 and the PDG valuei, 
412 ± 2(stat.) fs and 412.6 ± 2.8 fs respectively. 

However many systematic effects on lifetime cancel in the lifetime ratio and therefore in y. 

The systematic uncertainties in y were estimated using large signal simulated data samples 
that were modified by variations in the event selection criteria and suitable variations reflecting 
the current understanding of the detector, the uncertainties in background level and composition, 



Table 1: A summary of the systematic uncertainties in the y measurement. 





y Uncertainty (%) 


Systematic Uncertainty 


K-K+ 


7T 7T + 


Tracking 


0.2 


0.9 


Particle Identification 


0.2 


0.4 


Background & D* Fragmentation 


0.2 


0.6 


Alignment and Vertexing 


+0.2 
-0.1 


+0.3 
-0.1 


Monte Carlo Statistics 


+0.4 
-0.6 


+0.4 
-0.9 


Quadrature Sum 


+0.6 
-0.7 


+1.2 
-1.4 



Table 2: A summary of the y results. The first error reported is statistical; the second, systematic. 



Decay Mode 


y (%) 


K-K+ 

7T~7r + 


1.5 ±1.318;? 
1.0 ±1.7!^ 


average 


1.4 ±1.0 



in beam spot position and size. The y estimation in simulated data shows no bias within 
statistical errors. Systematic checks of the SVT internal alignment have been performed using 
e + e~ — >77~ >4 prongs events which provide a high statistics, zero lifetime control sample. All 
systematic uncertainties in the y parameter are summarized in Table 1. Major contributions 
to tracking uncertainties are derived by transverse momentum {pt) scale and tt s resolution 
uncertainties and, particularly for the tttt mode, by introducing a q/pT costant shift. Particle 
identification uncertainties are mainly associated to DIRC Cerenkov angle resolution (especially 
for the 7T7T mode) . Background uncertainties are introduced for the tttt mode by the level variation 
of the zero lifetime background and the uncertainty in the portion assigned to the background 
at lower mass values of the signal region. The major contributions to vertexing/alignment 
uncertainties are associated to a constant/partially additional smearing in the y position of the 
interaction point (a variation of beam spot size). However the largest source of uncertainty is 
from Monte Carlo statistic and will be certainly reduced in the future. 

6 Results and near term prospects 

The measured values of y for the two Cabibbo-suppressed decay modes separately are presented 
in Table 2 together with their average to be compared with the y limits of similar precision 
provided by BELLE0 and FOCUsi collaborations, (-0.5 ± 1.0±g;|)% and (3.42 ± 1.39 ± 0.74)% 
respectively. Our result is consistent with zero but suggests a positive value not incompatible 
with the FOCUS result. This interesting but by no means conclusive result induces us to perform 
in the future a new measurement over a rather larger data sample. Indeed the statistical error 
is determined by the amount of the Cabibbo-suppressed decay events selected. 

Limits on x = AM/T and y derived by the time-dependent analysis of D° — » K + ir~ wrong- 
sign decays are going to be released in the near future by BaBar. 
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